Reconnection and acoustic emission of quantized vortices in superfiuid 
by the numerical analysis of the Gross-Pitaevskii equation 
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We study numerically the reconnection of quantized vortices and the concurrent acoustic emission 
by the analysis of the Gross-Pitaevskii equation. Two quantized vortices reconnect following the 
process similar to classical vortices; they approach, twist themselves locally so that they become 
anti-parallel at the closest place, reconnect and leave separately. The investigation of the motion of 
the singular lines where the amplitude of the wave function vanishes in the vortex cores confirms that 
they follow the above scenario by reconnecting at a point. This reconnection is not contradictory to 
the Kelvin's circulation theorem, because the potential of the superflow field becomes undefined at 
the reconnection point. When the locally anti-parallel part of the vortices becomes closer than the 
healing length, it moves with the velocity comparable to the sound velocity, emits the sound waves 
and leads to the pair annihilation or reconnection; this phenomena is concerned with the Cherenkov 
resonance. The vortices are broken up to smaller vortex loops through a series of reconnection, 
eventually disappearing with the acoustic emission. This may correspond to the final stage of 
the vortex cascade process proposed by Feynman. The change in energy components, such as the 
quantum, the compressible and incompressible kinetic energy is analyzed for each dynamics. The 
propagation of the sound waves not only appears in the profile of the amplitude of the wave function 
but also affects the field of its phase, transforming the quantum energy due to the vortex cores to 
the kinetic energy of the phase field. 
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I. INTRODUCTION 

The numerical simulations of the quantized vortex dy- 
namics in superfiuid are classified into two methodstl. 
The first is the vortex filament formulation and the sec- 
ond is the numerical analysis of the Gross-Pitaevskii 
equatiou/GPE). The former simulation pioneered by 
Schwarza is useful for superfiuid 4 He, since the core of 
quantized vortex is very thin, i.e. of the order of atomic 
size. In this formulation, the motion of vortices is caused 
by the local and the nonlocal induced velocity fields sub- 
ject to the Biot-Savart law. When two vortices approach 
each other within a critical distance, they are assumed 
to reconnect. The results of this formulation agree excel- 
lently with the experimental results of the vortex tangle 
in superfiuid turbulence. However, the significant phe- 
nomena concerned with the core structure, such as vor- 
tex nucleation, annihilation and reconnection, cannot be 
described by this formulation. 

On the other hand, the analysis of the GPE, which 
is discussed in this paper, can investigate and analyze 
the superfiuid dynamics including the motion of vortex 
cores. Koplik and Levine solved numerically the GPE 
and showed two closed vortices actually reconnected even 
in an inviscid fluidtj. However, they did not show the de- 
tail of the reconnection process and the concurrent acous- 
tic emission, i.e. the emission of the condensate density 
waves. This is a significant phenomenon in the vortex 
dynamics. 



As well-known, the reconnection of vortices and the 
concurrent acoustic emission occur in a classical viscous 
fluiaD. In this case, the viscosity plays an important role 
in the reconnection. The deduction from this phenom- 
ena is not applicable to the inviscid superfiuid. The im- 
portant difference between classical and quantized vortex 
cores is the density profile. The quantized vortex has the 
core where the superfiuid density vanishes. On the other 
hand, the vortex core in an ordinary classical fluid has 
the finite fluid density. It is the core structure that causes 
the reconnection of quantized vortices even in an inviscid 
fluid. 

The reconnection and the concurrent acoustic emission 
in superfiuid are significantly related with the "eddy vis- 
cosity" which is a long-standing problem in superfiuid 
physics. The superfiuid turbulent stateu in a capillary 
flow induces excess temperature and pressure differences 
between both ends of the capillary. It is widely recog- 
nized that the temperature difference arises from the mu- 
tual friction due to the scattering of the excitations by 
the quantized vortex. On the other hand, the pressure 
difference is described phenomenologically by the eddy 
viscosity. The eddy viscosity works for superfiuid and 
reduces its total momentum, but its origin has not been 
necessarily revealed. The reconnection and the concur- 
rent acoustic emission reduce the energy of the quantized 
vortex, and then increase the kinetic energy of the super- 
fluid including the sound waves. If the kinetic energy de- 
cays, the total energy and momentum of the superfiuid 
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are reduced. Therefore, the reconnection and the con- 
current acoustic emission may be the origin of the eddy 
viscosity. _ 

Recently Leadbeater et al.a found that a sound wave 
is emitted when two vortex rings reconnect. They cal- 
culate the energy of the sound wave by estimating the 
reduction of the vortex line length. Our works investi- 
gate the change of energy components when two vortices 
reconnect, following the method of Nore et aEi. 

In this paper, we show the detailed process of the re- 
connection of quantized vortices by the numerical anal- 
ysis of the GPE, and the concurrent acoustic emission. 
Section [n] describes the basic equations and the numer- 
ical method. In Sec. Ill, we show the acoustic emis- 



sion when two anti-parallel vortices approach and dis- 
appear. This acoustic emission could be related with 
the Cherenkov resonances. The obtained solutions are 
concerned with the-stationary soli tary waves studied by 
Jones and RobertsO. In Sec. IV, to show the detailed 



process of the reconnection, we analyze the motion of 
two vortices which are initially placed at a right angle. 
Section [V| describes the sequence of reconnections and 
concurrent acoustic emissions. This process resembles 
Feynman's cascade processEj. Section VI describes the 



change in energy components. Considering the quantized 
vortex reconnection, we need to discuss whether it con- 
tradicts the Kelvin's circulation theorem that states the 
conservation law of the ci rcula tion in a n idea l fluid, which 
will be described in Sec. VII. Section VIII is devoted to 



conclusions and discussions. 



II. BASIC EQUATIONS AND NUMERICAL 
METHOD 

The GPE is 



dt 



2m 



V 2 ^{x) + g\^{x)\ 2 ^{x), 



(1) 



where ty(x) is the macroscopic wave function, and the 
chemical potential /i is renormalized by the global gauge 
transformation ^(x) = ^>'(x)exp(—i/j,t/fi). The coeffi- 
cient g of the interaction is related with the s-wave scat- 
tering length a as g — 4Trh 2 a/m. The GPE conserves the 
number of particles and the total energy given by 



N 



d 3 x\Hf(x) 



d 3 x \j-\V*(x)Y 



(2) 



(3) 



respectively. This wave function is expressed in terms of 
a density and a phase by Madelung's transformation 



= \J p(x) exp (iO{x)) . 



(4) 



The superffuid velocity is given by v(x) — TiV9/m, where 
m is the mass of bose particles. 

The total energy i? to t can be decomposed as the 
following!]: 



E u 



E\,\ 



kin 



E kin = / d 3 x— (Vp( x ) v ( x )) > 



E int = / d d xg(p(x)Y 



E *= MXz: vVp(* 



(5) 
(6) 

(7) 
(8) 



The kinetic energy E\a n is related with the velocity field, 
and £i n t is the internal energy of the fluid. The quantum 
energy E q comes from the gradient of the condensate, so 
that it is large when a vortex exists, as described later. 
Furthermore, in order to estimate the compressibility ef- 
fects, we decompose ^fpv into ^fpv = (^/pvY + (y/pv) c 
with V • {yfpv) 1 — 0. The corresponding components 
are named ££ in and E^ in , satisfying the relation Eki a = 
E^ in + E^ in . The compressible kinetic energy E^ in is di- 
rectly concerned with the acoustic emission. 

The numerical analysis is made for the normalized 
GPE. Equation (Q) is reduced to 



4[=-v 2 /+i/i 2 /, 



(9) 



where / = l 3 / 2 ^, x x/l, t -> t/{2ml 2 /K), I = 8vra, 
/ = an d the healing length £ = l/^/87ran . Small 

perturbations around the uniform stationary solution 
yield the sound wave of velocity c = \[2. A straight vor- 
tex is given by the stationary solution where |/| increases 
from zero at the central core over the healing length to 
unity at infinity. This two dimensional solution is ex- 
pressed as /o = l/MIe 1 ^ in the polar coordinates (r, <p). 
An initial vortex configuration is determined by its core 
location and vorticity direction. For example, a vortex 
with a core at (xotVo) in the direction of z axis gives 
/(#, y, z) = fo(x — Xq, y — yo). The initial configuration 
with several vortices is given by multiplying each of the 
single vortex. 

To numerically; integrate the GPE, we use a split step 
Fourier methodi3. There are several advantages of this 
method: (i) each time step is split into two segments, 
the first of which integrates the non-linear term in real 
space, and the second integrates the Laplacian operator 
in Fourier space, (ii) this method is second order accurate 
in the time interval of numerical analysis and all order in 
the space interval of that. 

This method should take careful account of the place- 
ment of vortices. First, we use the periodic bound- 
ary condition, because it allows us to use Fast-Fourier- 
transformation, which saves computational times. Sec- 
ondly, for example, a single vortex configuration in x- 
y plane is shown in Fig. [I]. Although we need only 
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[0, L] x [0, L] region, the periodic box [-L, L] x [-L, L] 
is used in the actual numerical calculation so that we 
can prevent the phase discontinuity of the vortex. If 
[0, L] x [0, L] may be adopted as the periodic box, the 
phase discontinuity occurs at the boundary between pe- 
riodic boxes. To prevent this discontinuity, the pe- 
riodic box must be [—L,L] x [—L,L], and we would 
add anti-vortices at (— xo,yo) and (xq,— yo), a vortex at 
(— Xo, —yo)- This procedure makes the phase well-defined 
everywhere. 



• 

(-xo,Vo) 


O 

Oo,yo) 




O 

(-xo,-yo) 


• 

(x ,-y ) 



-L 



-L L „ 

FIG. 1. Configuration of a single vortex on the [—L,L] 

periodic box. When we place a vortex at (xo,J/o), to prevent 

the phase discontinuity, we add a vortex at (—so, — yo) and 

anti-vortices at (— Xo,yo) and (xo, —yo)- 

In the following figures, we usually show only [0, L] 3 
part for three-dimensional calculation and [0, L} 2 part for 
two-dimensional one. The conservation of the normaliza- 
tion and the total energy is confirmed numerically. These 
quantities are conserved in the order of 10 -6 in all simu- 
lations. 



III. ACOUSTIC EMISSION IN THE 
ANNIHILATION OF TWO ANTI-PARALLEL 
VORTICES 

Consider the initial configuration where two anti- 
parallel vortices are placed as shown in Fig. |^(a). This 
is a simple and typical configuration, whose develop- 



ment includes the Che 
tices studied by Ivonin 
studied by Jones et al 



;nkov resonance of quantized vor- 
and the stationary solitary waves 
. The former is deeply concerned 
with the acoustic emission when two anti-parallel vor- 
tices approach and disappear. The latter shows the an- 
alytical solutions of the GPE. Jones et, al. investigated 
the axisymmetric stationary solution of the GPE to ob- 
tain a continuous family consisting of two branches in the 
momentum- energy plane. One is a vortex ring(an anti- 
parallel vortex pair), and the other is a rarefaction pulse 
without vorticity. 



As described in Appendix, two anti-parallel vortices 
approach each other under the periodic boundary condi- 
tion. When the distance between them is reduced within 
a critical distance 2£, the velocity of two vortices be- 
comes comparable to the sound velocity, then the lo- 
cal Chetenkov resonance starts and the sound waves are 
emitted!! The energy of the vortices is partly lost to 
the sound emission, the distance between them being re- 
duced. These processes continue until two vortices dis- 
appear. The sound wave emitted at the moment of the 
annihilation has the largest amplitude in the process of 
the Cherenkov resonance. 

In Fig. 0, the initial distance between two anti-parallel 
vortices, ALq = 1.5625£, is less than 2£, so that they 
approach each other(Fig. ||(b)). Their singular lines of 
|/|=0 overlap for a very short time and disappear. 




FIG. 2. Annihilation of two anti-parallel vortices, t = 0(a) 
and t — 0.242(b). The surface represents the contours of 
|/| =0.06. 

We consider the change of the condensate density |/| 
in the process of Fig. ^|. Note that it is necessary to 
separate carefully the effect of acoustics from that due to 
vortex motion. In order to reveal both the vortex motion 
and the acoustic emission, we show the change of |/| in 
Fig. H Initially two very close depressions represent the 
singular cores of the vortices ((a)). They become close 
and overlap for a very short time and disappear ((b)). 
After their annihilation at t ~ 0.2, the vestigial of two 
vortex cores becomes flattened gradually((c)-(f)). The 
sound propagations are dimly shown in Fig. 0; in order to 
show them more clearly, Fig. ^| shows the sectional profile 
of the density |/|. After the annihilation, some wakes are 
found to propagate outside over the global profile; this is 
nothing but a sound wave. Its propagation velocity is of 
the same order as c = V2- 

The motion of the vestigial of two vortex cores is re- 
lated with the stationary solitary wave studied by Jones 
et an. We investigate the contour plots of the condensate 
density. These plots resemble those obtained by Jones et 
al, although they studied the stationary solution. Our 
dynamical solution follows approximately the vortex so- 
lution, makes two vortex cores disappear, then follows the 
rarefaction pulse solution. The vestigial of the vortices 
correspond to the rarefaction pulse. The above sound 
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waves propagate over the profile of the rarefaction pulse. 




FIG. 3. Profile of |/| when two anti-parallel vortices col- 
lide. The vertical axis refers to |/|, and the others represents 
the space coordinates, t = 0(a), t = 0.242(b), t = 0.324(c), 
t = 0.61(d), t = 0.814(e), t = 1.14(f). The depressions 
disappear(a)-(b), and the core structure are flattened grad- 
ually^)-^). The dimly seeable wakes are sound waves(f). 




FIG. 4. Sectional profile of |/|. These plots show onl 
0.895 < |/| < 0.915. Coordinates are the same as Fig 
t = 0(a), t = 1.14(b), t = 1.18(c), t = 1.222(d). Some propa- 
gated wakes are sound waves. 



The change in each energy component is shown in Fig 
M This figure does not show the internal energy com- 
ponent -Eint, because it is much larger than the other 
components, E- lnt /-Etotai ~ 1- The compressible and in- 
compressible energy increases after two vortices vanish at 
t ~ 0.2. The increase in the compressible energy is rea- 
sonable. On the other hand, the increase in the incom- 
pressible energy may remain controversial. The quantum 
energy of the vortex cores is transferred irreversibly into 
the kinetic energy as described later. 
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FIG. 5. Change in each energy component for the process 
of Fig. []. The thick broken line , the thin broken line and 
the dotted line shoe E q /Etot, 



-Ekin/ El 



lot and E^ in /E t , 



The sound wave affects not only the density |/| but 
also the field of the phase 9. The vector field v s = V6* in 
x-y planes are shown in Fig. @. The symbols and (g> 
denote the singular lines with \ f\ = of a vortex and an 
anti-vortex, respectively. Two vortices become close((a)- 
(c)) and disappeared)). Before their annihilation , they 
become close together with the circulative field around 
itself. Just before the annihilation, \v s \ between two vor- 
tices grows very much but the density |/| there is reduced. 
When their singular lines overlap, v s between them dis- 
appears, while the circulative field around them still re- 
mains((d)). Since the propagation speed of the event, 
i.e. "the annihilation of the singular lines", is of the 
same order as c, v s far apart from the annihilation point 
keeps the circulative flow for a certain time((d)-(f)). Af- 
ter this event information comes up, the circulative flow 
disappears((g)-(h)). The arc shown in Fig. ^| indicates 
the crest of the sound wave emitted at the annihilation, 
which is estimated by the density propagation as shown 
in Fig. H and [|. The arc shows more clearly the disap- 
pearing processes of the circulative flow. As described 
before, the quantum energy of the vestigial of the sin- 
gular cores are reduced and transferred into the kinetic 
energy. The increase in the kinetic energy results from 
the large velocity field near the annihilation point shown 
in Fig.f ff)-(h). 

Figure fi] shows the change of the vector field from the 
initial state 
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Sv(x,y,t) 



\v(x,y,t) - v(x,y,Q)\ 
\v(x,y,0)\ 



(10) 



Ml 



These figures show the propagations of the change of the 
vector field, reflecting to the propagation of the sound 
waves as described above. The distance between two 
anti-parallel vortices oscillates by the acoustic emission 
of the Cherenkov resonance!, including the contribution 
from the vortex motion too. Although the circulative 
flow is slightly disarranged by the sound wave by the 
Cherenkov resonance for a moment, then it recovers the 
previous state((f)-(g)) mostly. This recovering mech- 
anism breaks down gradually as the amplitude of the 
sound wave increases, eventually the circulative flow dis- 
appears after the maximum sound wave emitted at the 
annihilation passes by. 

We investigate also the distribution of the incompress- 
ible kinetic energy density 



and the compressible kinetic energy density 

e kin( X ) = {(VP") ) 2 . 



(11) 



(12) 



The distribution of e£ in (a:) has small peaks on the sound 
waves and is also spread into the whole system. On the 
other hand, that of e£ in (a;) has a very large peak near 
the vestigial of the singular cores, which contributes to 
the increasing in the incompressible energy in Fig. ^. 
Eventually, because there is no dissipative mechanism, 
these energy components are mixed and oscillates after 
that. 
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FIG. 6. Vector field V0, t=0(a), t=0. 14(b), t=0. 192(c), 
t=0.316(d), t=0.406(e), t=0.61(f), t=0.814(g), t=1.222(h). 
The symbols and ® denote a vortex and an anti-vortex, 
respectively. The vector heads denote the direction of V#, 
and its length denotes |V#|. The symbol * shows the annihi- 
lation point, and the arc shows the crest of the sound wave 
emitted at the annihilation. For convenience of visualization, 
the vector length is stretched. 
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FIG. 7. Change 5v(x,y,t) of the velocity field from the 
initial state, t=0(a), t=0. 14(b), t=0.192(c), t=0.316(d), 
t=0.406(e), t=0.61(f), t=0.814(g), t=1.222(h). 



IV. DETAILED PROCESS OF RECONNECTION 

In this section, to show the detailed processes of the 
reconnection, we study the three-dimensional dynamics 
starting from the configuration where two vortices are 
placed at a right angle as shown in Fig .||(a). The grayed 
surfaces represents the contour of |/| = 0.1, which is more 
inside in thet-icore than that of |/| = \/0.3 calculated by 
Koplik et aB. Two vortices cause local twists in each 
other((b)) so that they become locally anti-parallel at 
the closest point((c)) and approach each other, then re- 
connected)) and leave separately((e)-(f)). This motion 
is the same as that found by the filament calculation!! 
The locally anti-parallel vortices which are closer than 
the critical distance become to follow the scenario of the 
Cherenkov resonance described in Sec. 



Ill: the anti- 



parallel part moves locally with the velocity comparable 
to the sound velocity, emits the sound wave and crosses at 
a point, corresponding to the two-dimensional annihila- 
tion. In order to investigate the switching of the singular 
core, we check the motion of the surface of |/| = 0.03 
too as shown in Fig. ||. Two lines in Fig. ^(a)-(d) show 
the singular line of |/| ~ 0. These lines reconnect at 
a point, i.e. the reconnecting point. Two vortices are 
found to reconnect at t ~ 1.45. The detailed processes of 
the reconnection are as follows: the singular cores of two 
vortices cross once, and keep the transient crossing for a 
very short time, path through each other, eventually leav- 
ing by changing their topology. The quantized vortices 



reconnect without going through the classical bridge-like 
reconnection processes. The quantized vortex reconnec- 
tion is not contrary to t he Kelvin's circulation theorem, 
as discussed in Sec. 



VII. 
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FIG. 8. Reconnection of two vortices, t = 0(a), 0.816(b), 
1.224(c), 2.448(d), 2.856(e), 3.06(f). The surface represents 
the contour of |/| = 0.1. Two vortices cause local twists 
in each other(b) so that they become anti-parallel(c) , then 
reconnect (d) and leave separately (e),(f). 
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(b) 



/ 



V 



(c) 



(d) 



FIG. 9. More detailed reconnection process of two vortices, 
t = 1.386(a,b), 1.428(c), 1.46(d). The surface represents the 
contour of |/| — 0.03. Two lines show the singular line by 
tracing the local minimum density less than |/| ~ 0.06. 

Figure [l0| shows the changes of the ratio of each energy 
component of the process of Fi g. |s| , which does not show 
Ei n t by the same reason in Sec. Ill . During the processes 
of the reconnection, E^ in and El in are reduced, and E q 
increases. If there would be the concurrent acoustic emis- 
sion, the compressible component should increase. The 
change of the energy components will be discussed in de- 
tail in Sec. VI. 
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V. SEQUENTIAL PROCESSES OF 
RECONNECTION AND CONCURRENT 
ACOUSTIC EMISSION 

Since the sound wave propagates over the density pro- 
file of the vortex and the vortex itself moves, we may find 
more easily the sound wave propagation in the region 
where the density of the condensate is nearly homoge- 
neous. Thus in order to find the sound wave propagation 
after the reconnection of 90° vortices under the periodic 
boundary condition, we shift the reconnecting point to- 
ward the corner of the periodic box in this section. The 
initial configuration is shown in Fig. [ll](a) . The vortices 
are placed at a right angle, one has the singular core along 
— z direction at (x, y) = (2.5, 1.9), other has that along y 
direction at (x, z) = (1.3, 1.9). In this configuration, we 
may show the density propagation enough far form the 
singular core. 

Figure 1 1 shows the dynamics of two very close vortices 
meeting at a right angle. These figures show that the vor- 
tices are reduced to the vortex rings through the recon- 
nections and disappear eventually. This processes is the 
cascade process proposed_by Feynmar£3 and revealed re- 
cently by Tsubota et. ali3. For example, the vortices are 
broken up to smaller loops by the reconnection((d)-(e)). 
A small one makes self-reconnection and disappears((f)- 
(g)), when the energy of vortices is transferred into the 
kinetic energy. This is just the final stage of the cas- 
cade process. At very low temperatures, these mecha- 
nism should play a crucial role in the decay of the vortex 
tangle. 

The changes in the ratio of each energy component are 
shown in Fig. [l2[ The figures below the horizontal axis 
indicate when the reconnections occur. 

The density propagation is visible enough far from the 
vortex core. The dimly seeable wakes are shown in Fig. 
|j~3| . In Fig. |lj which enlarges the sectional profiles of 
Fig. [13], the propagation of the wakes are shown prop- 
erly. From the analysis of the speed of the sound wave in 
Sec. Ill , we find that the sound wave © and probably 



comes from the reconnection @ and @, respectively. As 
described before, the amplitude of the sound wave de- 
pends on how much the quantum energy is transferred 
into the kinetic energy, in other words, the vortex line 
length is lost by the reconnection. The reconnection © 
is the point reconnection, while @ is the annihilation of 
a small vortex ring. Although the lost line length of each 
two reconnection is shorter than that in Sec. Ill, the 



sound wave propagations are obviously visible. This is 
because the density enough far from the singular cores 
is about homogeneous. For example in 5 < x < 10 and 
5<2/<10, 1/1-0.8. 



FIG. 10. Change in the energy component for the pro- 
cess of Fig. ^[ The thick broken line ,the thin broken line 
and the dotted line show respectively E q /E to t, El. in /E to t and 

-Skin/^tot. 
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FIG. 11. Reconnections of two vortices, t = 0.0(a,b), 
0.182(c), 0.58(d), 0.928(e), 1.326(f), 1.794(g), 2.132(h), 
3.498(1), 3.948(j), 5.172(k). The surface represents the con- 
tour of |/| = 0.06. The inside box of (a) indicates the region 
plotted in (b)-(k). 





FIG. 12. Changes in each energy component of Fig. [n]. 
The thick broken line ,the thin broken line and the dotted 
line show respectively E q /E to t, El in /E tot and El in /E to t- The 
numerics on the bottom show when the reconnection occurs. 




V 
- U.S2 
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FIG. 14. Sectional profile of |/|. The range of |/| is limited 
to 0.76 < |/| < 0.82. Coordinates are the same as Fig. [l^, 
t = 0.0(a), 1.938(b), 2.04(c), 2.142(d). Some propagating 
wakes are sound waves. 



VI. CHANGE IN THE ENERGY COMPONENTS 



-4? 




(c) 



(d) 



FIG. 13. Profiles of |/|. The vertical axis refers to |/|, 
and the others represents the space coordinates, , t = 0.0(a), 
1.938(b), 2.04(c), 2.142(d). The dimly seeable wakes are 
sound waves((b)-(d)). 



The comparison between Fig. § Fig. [h] and Fig. 
|l2| shows some important properties on how the energy 
components behave. The energy is generally transferred 
between these three components. In fact, the changes in 
the components after the reconnection are roughly clas- 
sified into two stages. In the first stage, E q increases, as 
shown in Fig. |o|, and the period < t < 1.2 of Fig. 
Then, the vortices reconnect at points and still survive 
after that. The second stage is characterized by the re- 
duction of E q in the period 0.2 < t of Fig. ||, and the 
period 1.2 < t of Fig. lj: the vortex cores rather collide 
and disappear, than reconnect at points. 

The behavior of the energy components in the first 
stage may be attributed at the following interaction be- 
tween the vortices and the sound waves. While two vor- 
tices reconnect at a point, as shown in Fig. ||, they stretch 
and fat near the reconnection point because of the inter- 
vortex interaction. The vortices near the reconnection 
point move with the velocity comparable to the sound 
velocity. The sound waves made by the reconnection, of 
course, expand outside, but the above factors make the 
sound waves interact continuously with the vortex cores 
even after the reconnection. Although we do not know 
the detail of the interaction between quantized vortices 
and sound waves, this continuous interaction must make 
it easy that the sound waves are absorbed by the vortex 
cores, which increases E q . 

On the other hand, the vortex cores fade out in the 
second stage. Thus the quantum energy decreases, and 
the kinetic energy increase owing to the acoustic emis- 
sion. However, Vp never vanishes even after the core 
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structure disappears, so that the energy components are 
mixed and oscillate. 



VII. KELVIN'S CIRCULATION THEOREM 

This section comments on the conservation of the cir- 
culation r = § ds ■ v, which is called the Kelvin's circu- 
lation theorem. This theorem states that the circulation 
round any closed material loop is invariant in an invis- 
cid(barotropic) fluid, and vortices are never born, die and 
reconnect. It seems that the theorem does not allow the 
above reconnection of quantized vortex. Then, we have 
to consider two problems: (i) is this theorem applicable to 
the superfluid dynamics derived from the GPE? (ii) what 
happens to the material loops when quantized vortices re- 
connect? Equation (||) yields the equation of motion for 
superflow: 



Dv 
Hi 



v 2 l/l 

l/l 



Vw, 



(13) 



where D/Dt = d/dt + v ■ grad. By differentiating T with 
respect to t, we obtain 



— - = <b ds ■ V \ w H 

Dt J V 2 



(14) 



Therefore T = const, if w is well-defined on the mate- 
rial loop along which T is defined((i)). When two vor- 
tices are about to reconnect, almost all material loops 
are carried away from the merger region by the flow(Fig. 
|l5|(a))((ii)). The loop just going through the reconnect- 
ing point breaks this theorem because w is not well- 
defined owing to |/| = 0, thus leading to the reconnec- 
tion(Fig. [l5](b)). This is the peculiar situation in super- 
fluid. Thus the reconnection of quantized vortices is not 
contradictory to the Kelvin's circulation theorem. 



Reconnecting point 




FIG. 15. Motion of the material loops when two vortices re- 
connect. (a):most loops are carried away from the merger re- 
gion by the flow. (b):the Kelvin's circulation theorem breaks 
on the loops going through the reconnecting point. 



VIII. CONCLUSIONS AND DISCUSSIONS 

We studied numerically the reconnection of quantized 
vortices and the concurrent acoustic emission by the anal- 
ysis of the GPE. The energy components are analyzed for 
every dynamics. 

When two anti-parallel vortices are close within a crit- 
ical distance which is twice the healing length, they move 
with the velocity comparable to the sound velocity, emit 
the sound waves and disappear; this phenomena is closely 
concerned with the Cherenkov resonance. The propaga- 
tion of the sound waves not only appears in the density 
profile of the wave function but also affects the phase 
field; they destroy the circulative velocity field that be- 
longed to the vortices originally. 

Two straight vortices meeting initially at a right an- 
gle follow the scenario similar to the vortex filaments, 
thus reconnecting. The singular lines of |/| = cross 
once at a point, switch their topology instantly and leave 
separately. This reconnection is not contradictory to the 
Kelvin's circulation theorem, because the potential of the 
superflow field becomes undefined at the reconnection 
point. 

Two close straight vortices meeting at a right angle, 
under the periodic boundary condition, follow a series of 
reconnection, breaking up to smaller vortex loops. Even- 
tually these loops disappear through self-reconnections 
with the acoustic emission. The process transforms irre- 
versibly the quantum energy of the vortex cores to the 
kinetic energy. This may represent the final stage of the 
vortex cascade process. 

The change in energy components depends on how the 
vortices reconnect. When vortices survive still after re- 
connections, the emitted sound waves interact continu- 
ously with the vortex cores. Then sound waves are ab- 
sorbed by the vortex cores, which increases the quantum 
energy. On the other hand, when vortex loops disappear 
through reconnections, the quantum energy is transferred 
into the kinetic energy with the acoustic emission. 

Throughout this work, it remains unresolved how 
sound waves interact with quantized vortices; they may 
stretch-, and fat the vortex cores or excite the Kelvin 
waveal3. More detailed studies of the interaction are be- 
ing carried out. 
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APPENDIX: MOTION OF ANTI-PARALLEL 
VORTICES BY THE FILAMENT FORMULATION 

This appendix describes the motion of the vortices ini- 
tially placed as shown in Fig. [l6| by the vortex filament 
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formulationO. The velocity field generated by a vortex 
whose singular core is at the origin is given by 



K 

2nr ' 



(Al) 



where k is the quantized circulation. The motion of equa- 
tion of the vortex 1 is given by 



s = Vx = v 2 i +v 31 +v 41 , 



(A2) 



where s = (x,y), and V21, U31 and v^i are the velocity 
field at the position of vortex 1 generated by the vortex 
2, 3 and 4, respectively. Each velocity are given by 



V21 



2tt 



0, 



V31 



K 1 

(sin ( 

2tt 2r y 



W41 



K 

27 



1 

2x 

COS 1 

1 



2y 







(A3) 
(A4) 
(A5) 



Therefore, the dimensionless form of the equations of mo- 
tion are 



y 



y 



1 

X 



y 



(A6) 



where t ■ k/2tt — > t. These equations can be rewritten in 
term of the polar coordinates: 



r = — cot 20, 
r 



1 

= ->0. 

r 



Therefore 



dv 

— = -2rcot26>. 
d9 



(A7) 



(A8) 



For example the vortex 1 is initially placed at (r, 0) 
(r ,7r/4), 



r 



sin 26. 



(A9) 



This equation can be rewritten in term of the Cartesian 
coordinates: 



(A10) 



This solution means that the vortex pair above x axis 
moves upward and approaches each other down to a cer- 
tain distance determined by the initial configuration as 
shown in Fig. [l^ (b), i.e. AL = tq. 

If AL becomes comparable to the healing length, the 
vortex filament formulation becomes not available, lead- 
ing to the pair annihilation as described in the text. 



I:- 




(a) 



4 



(b) 



FIG. 16. (a) The initial configuration of SecjlEj. The © 
and ® symbols indicate a vortex and an anti-vortex, (b) The 
trajectories of each vortex, the vortex pair above x axis moves 
upward and approaches each other down to a certain distance. 
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